Capillary condensation for fluids in spherical cavities 
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The capillary condensation for fluids into spherical nano- 
cavities is analyzed within the frame of two theoretical ap- 
proaches. One description is based on a widely used simplified 
version of the droplet model formulated for studying atomic 
nuclei. The other, is a more elaborated calculation performed 
by applying a density functional theory. The agreement be- 
tween both models is examined and it is shown that a small 
correction to the simple fluid model improves the predictions. 
A connection to results previously obtained for planar slits 
and cylindrical pores is done. 
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PACS numbers: 61.20.-p, 64.70.-p, 67.70. +n 



I. INTRODUCTION 

Since more than a decade a so called simple fluid model 
(SFM) is being applied for analyzing wetting properties 
of classical and quantum fluids adsorbed onto planar sub- 
strates [1,2]. The chief idea of this description is to write 
down the free energy of the adsorbate as a sum of domi- 
nant volume, surface, and substrate terms. As shown in 
Refs. [3] and [4], the use of this SFM for studying capil- 
lary condensation in the case of adsorption between two 
parallel planar walls (i.e. in a slit geometry) yields a rich 
pattern of phase transitions in agreement with more elab- 
orated calculations carried out with density functional 
(DF) approaches. More recently, the SFM has been ex- 
tended to investigate systems with cylindrical symmetry 
[5-9] . The authors of Ref . [7] have found that for a rather 
compact regular array of infinitely long, solid, and par- 
allel cylinders, besides the well known film and capillary 
condensation (CC) phases, an additional necking config- 
uration may appear. On the other hand, in Ref. [8] we 
have applied this SFM to describe CC transitions of no- 
ble gases confined into single cylindrical pores of alkali 
metals finding a nearly universal behavior. A detailed 
comparison of SFM and DF results performed in Ref. [9] 
indicates a fair agreement for the main features of cylin- 
drical pores. The phenomenon of hysteresis has been also 
addressed in Ref. [8] , it was found a good correspondence 
with results reported in the pioneering work of Cole and 
Saam [10]. 

Currently there is an increasing interest for investi- 
gating the storage capacity of small pores with differ- 
ent geometries (e.g., nanotubes and rather spherical zeo- 



lites, see Ref. [7] and references quoted therein). Accord- 
ing to this trend, the study of adsorption in spherical 
nano-cavities becomes relevant. Therefore the aim of the 
present work is to analyze phase transitions for classical 
and quantum fluids of noble gases confined into alkali 
metal cavities. We selected these systems because the 
involved interactions are known. Its behavior has been 
obtained from the SFM outlined above. In order to check 
the validity of such a rather crude model for this geome- 
try, the adsorption of 4 He is also calculated by applying 
a DF. The latter approach has been recently successfully 
utilized for investigating properties of 4 He adsorbed on 
very small spherical balls (with radii of a few A) of dif- 
ferent materials [11,12]. However, for the sake of com- 
pleteness, we should mention that Monte Carlo [13-15] 
as well as DF [15,16] calculations have been performed 
for classical gases adsorbed on planar alkali metals. 

The paper is organized in the following way. In Sec. 
II we describe the evaluation of the adsorption potential 
and outline both the SFM and the DF approaches for 
capillary condensation. The numerical results and an 
improvement of the bare SFM are presented in Sec. Ill 
together with a discussion of the hysteresis cycles and a 
summary of the whole picture exhibited by the systems. 
Section IV is devoted to final remarks. 



II. THE MODEL 

The properties of a fluid adsorbed by a solid substrate 
may be studied by analyzing the grand free energy 



Cl[N] =F[N}- [iN , 



(2.1) 



where F is the Helmholtz free energy, \x the chemical 
potential, and N the number of particles of the adsorbate. 
Quantity F[N] contains the energy due to the interaction 
between fluid atoms as well as the energy provided by the 
confining potential. Let us first outline our approach for 
the description of the fluid-substrate interaction. 



A. Adsorption potential 

In order to construct the physisorption potential inside 
a spherical cavity we assume that an adsorbed atom lo- 
cated at r, measured from the pore center, interacts with 
an elementary substrate volume at fl via an isotropic 
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(12,6) Lennard- Jones (LJ) pair potential with standard 
parameters Elj and ctlj. Next, we suppose that sub- 
strate atoms are uniformly distributed over the volume 
outside the cavity, this continuous substrate approxima- 
tion leads to a spherically averaged potential. Under 
these conditions, when the substrate extends from r = R p 
up to r = oo, the total potential may be expressed as 
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Here, p su b is the averaged density of the substrate in 
number of atoms per unit volume. Of course, each ap- 
proximation performed along this procedure introduces 
an error. However, we expect that the resulting poten- 
tial would give a fair description of the main features of 
examined systems. Notice that this procedure has been 
previously adopted in Refs. [17,18] for getting the plain 
(9,3) potential in the case of planar surfaces. 

After some straightforward algebra the integral over 9 
may be cast into the form 
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(2.3) 



where M n {v) stands for the integral 
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with v = r/ri. The integration over rl yields the follow- 
ing expression for the adsorption potential 
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(2.5) 



This formula may be compared with expressions derived 
in Ref. [19]. On the other hand, in the limit of a very 
large cavity (Rp — ► oo) an expansion in terms of z = 
R p — r << R p reduces Eq. (2.5) to the (9-3) potential 
corresponding to a planar substrate (cf. Eq. (4) in [18]) 
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(2.6) 



For all the inert gases, the substrate potentials pro- 
duced by pores in different alkali metals were calculated 
using values of Elj and <jlj determined by adjusting the 
pair potentials of Patil [20] according to a procedure out- 
lined by Ancilotto et al. [21]. Density p su b was evaluated 
by taking into account that solid alkali metals are bec 
crystals with a lattice constants a listed in Table 3 of 
[22]. All the adopted parameters are quoted in Tables I 
and II. 



B. Simple fluid model 

An expression for the Helmholtz free energy per par- 
ticle for the CC phase, /cc> m ay be obtained from the 
formalism developed in [5,6]. Assuming that the free en- 
ergy of the fluid is given by the main contributions of the 
volume, surface, and substrate terms; then, if the fluid 
fills a pore forming a sharp sphere of radius Rq with den- 
sity equal to the bulk equilibrium value po one gets 
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(2.8) 



Here /oo is the asymptotic value and it coincides with the 
chemical potential po at saturated vapor pressure (SVP). 
This expression for the free energy has been recently uti- 
lized for studying adsorption into cylindrical pores [6,8]. 

In all the cases, for a spherical cavity with a certain ra- 
dius R p the potential U su ^,(r) given by Eq. (2.5) exhibits 
a minimum with depth Db, v at a distance r m (R p ) from 
the wall. It can be analytically shown that for a given 
fluid-substrate combination, in the regime R p > <jlj, the 
depth Dr p decreases for increasing size of the spherical 
cavity and attains an asymptotic value Dqo according to 
the expression 
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Here r m is the asymptotic position of the minimum 
[r m {R p oo)] which is given by 



1/6 



(2.12) 



The behavior of D Rp /D^ as a function of 1/ (R* - 1/2) 
is similar to that previously obtained for the cylindrical 
geometry (see Figs. 1 and 2 in [8] and Fig. 1 in [30] 
). In order to facilitate any comparison the values of 
r m and are listed in Table II. Since the properties 
found for the adsorption into cavities in Rb and K do not 
differ much from those for pores in Cs, the correspond- 
ing values are not included in this table. Note that for 
a given fluid-substrate system, the present r m is larger 
than the corresponding parameter for a slit geometry z m 
determined from a (9,3) potential in Ref. [3] (cf. Table 
I therein). This feature is simply due to different ap- 
proaches used to fix the potential parameters. 



1. Capillary condensed droplet phase 

Upon taking into account Eqs. (2.7) and (2.8), and 
defining the sharp sphere radius i?o as the effective pore 
radius R n — R p —r m , the grand thermodynamic potential 
for the CC phase becomes 



^cc = F cc -pN 
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Starting from this equation, it is possible to get the mod- 
ified Kelvin equation (KE) for CC. By setting Occ = 0, 
one arrives at 

<J s i - ^ (At - Mo) Po (R P - r m ) = , (2.14) 

where a s i is the liquid-substrate interfacial tension de- 
fined as the sum 
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Next, by using the Gibbs-Duhem [31] relation for the 
fluid and neglecting any compression of the system one 
gets a link to the pressure P 

Po (p-po) = P-Po, (2.16) 

which leads to the usual KE 

3 a sl - (P - P ) {Rp - r m ) = . (2.17) 

This equation expresses the pressure reduction for con- 
densation in terms of the effective radius of the condensed 



fluid {R p — r m ). The purpose of this work is to generalize 
the KE by taking into account both: (i) the explicit de- 
pendence of f^cc on the substrate potential and (ii) the 
role of film formation on the adsorption behavior. 

For the sake of generality, we shall discuss the predic- 
tions in terms of the dimensionless grand potential 
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It is convenient to cast the integrated adsorption poten- 
tial per unit area Ijj in the following reduced form 
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Here, the strength parameter is the reduced asymptotic 
well depth 
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and g{R* — 1} is the dimensionless integral 



g{£} = 4 / t' 2 [-^.ub(r ro SO/D00] , (2-22) 
s Jo 



with 



C = r / r n 



(2.23) 



Consequently, one has also to introduce the reduced dif- 
ference of chemical potentials 



A = (p - p) r m po/cri v 



(2.24) 



These definitions of D* and A are formally equal to those 
adopted in previous works [3,8], provided that for planar 
parallel walls instead of r m one uses z m . 

Finally, one may write the reduced version of the grand 
potential given by Eq. (2.19) in the following way 



n 



cc = l--D*g{R*-l} + -A{R* -1) . (2.25) 



This expression allows a study of adsorption in terms of 
D* . Before going ahead, let us remind that a fundamen- 
tal property of the integral g{£} introduced by Eq. (4) 
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in [3] is its independence of the fluid-substrate combina- 
tion. This feature is also exhibited by the integral g{£} 
calculated according to Eq. (2.22) 
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For £ = i?* — 1, in the limiting case R* — > oo one gets 

If 



lim o{.R* - f } = 



(2.27) 



This result is equal to that obtained for the slit geometry 
[3]. 



2. Shell-film phase 

Under the same assumptions adopted for writing Eqs. 
(2.7) and (2.13), the grand thermodynamic potential for 
a spherical film of thickness £, i.e., for the shell-film (SF) 
phase, may be expressed as 

«sf = F SF -pN 

= 4:irai v [{R p - r m f + (R p - r m - if] 
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Note that the film grows from r = r m towards the center. 
The reduced version of this grand free energy reads 
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where x is the dimensionless thickness 

X = 1 1 Tra , 
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(2.29) 



(2.30) 



which is related with the dimensionless inner radius y 
measured from the center of the sphere 



y = {Rp-r m - l)/r m = R* - 1 - x . 



(2.31) 



In this case, in order to find the stable configuration at 
fixed R* (i.e., R p = const.), one must determine the value 
of x which provides the minimum fig F . 



C. Density functional theory 

In the DF approach the ground-state energy of an in- 
teracting ./V-body system of 4 He atoms, confined by an 
adsorbate-substrate potential U su b(r), may be written as 

#gs / dr v /^(r)V 2 v /^(r)+ J dr p(r) e sc (r) 

drp(v)U suh (r) , (2.32) 



where p(r) is the one-body density and e sc (r) the self- 
correlation energy per particle. The density profile p(r) 
is determined from the Euler-Lagrange (EL) equation de- 
rived from the condition 
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(2.33) 



In the case of a spherical symmetry the variation of Eq. 
(2.33) leads to the following Hartree like equation for the 
square root of the one-body helium density 
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V H (r) + U suh (r) 



VPM=MV / PW, (2-34) 



which also determines \x. Here Vh(t) is a Hartree mean- 
field potential given by the first functional derivative of 
the total correlation energy E sc [p] 
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The expression for the spherically symmetric Vn(r) de- 
rived in the case of the Orsay-Paris nonlocal DF (OP- 
NLDF) proposed in [32] is given in the Appendix of [33] . 
Equation (2.34) was solved at a fixed number of helium 
atoms 
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III. GENERAL RESULTS - PHASE DIAGRAMS 

A. Thresholds for CC and SF phases at SVP 

The case of SVP (A = 0) leads to a very simple cri- 
terion for the occurrence of CC. The expression for the 
transition line (Qqq — 0) which separates the behavior 
into two regimes, capillary condensation (Qq C < 0) or 
empty (n* cc > 0), is 



£>A=o(CC) = 



g{R* - 1} 



(3.1) 
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It coincides with Eq. (7) of [3] and with Eq. (3.1) of [8]. 
Its dependence on R* is displayed in Fig. 1. Therefore, 
one can state that the curve given by Eq. (3.1) provides a 
universal relation for the critical values of the parameters 
R* and D* A=Q . On the other hand, it should be noticed 
that at SVP the transition from empty (E) to SF would 
occur for 



^A=o(SF) 



2(i?*-l) 2 + 2(i?*-l-a;) 



(R* - l) 2 g{R* - 1} - (R* - 1 - x) 2 g{R* - 1 



x}- 
(3.2) 



It is clear that, for < x < R* — 1, Z?2 i=0 (SF) is always 
larger than D^ =0 (CC). So, for increasing D* at SVP, 
before Og F becomes zero Qq C is already negative favoring 
the CC phase against the SF one. 

Furthermore, the threshold condition for CC to occur 
at SVP for large R* changes very little for increasing R* 
and for very broad pores attains the asymptotic value 



D* A=0 (CC;R* -► oo) = lim 

_R*— »oo 



g{R* - 1} 
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n 



t 2.9. 
(3.3) 



A glance at Table II indicates that noble gases heavier 
than 4 He adsorbed in cavities of Cs would not form in the 
present approach a stable CC phase at the triple-point 
temperature T = T t , because the corresponding values of 
D* are smaller than 2.9. 



B. Thresholds for CC and SF phases below SVP 

The general problem of behavior below SVP is more 
complicated. The presence of the A term in Eq. (2.25) 
leads to the following threshold value for CC 



D*(CC) = D£ =0 (CC) 



l + -A(iT-l) 



(3.4) 



However, one must now examine the possibility that a 
film configuration has a lower free energy than that of 
the CC phase. So, in the case A > 0, it is necessary to 
evaluate the minimum of f2g F as a function of the variable 
x and to compare the result with f^c- 

The reduced depth for E — > SF transition at A > 
may be obtained from Eq. (2.29), written in terms of y 
it becomes 



D*(SF) = £>1 =0 (SF) 
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(3.5) 



The difference between the grand free energy of the 
SF phase and that of the CC case written in terms of the 
dimcnsionless inner radius y is 



cc 



R* 



±D*g{y}-±Ay 



(3.6) 



Here, in the second parenthesis, the first term represents 
the extra surface energy of the film, the second provides 
the interaction between the solid substrate and the ad- 
sorbed atoms which fill the gap when the SF to CC transi- 
tion occurs; while the third stands for the free energy cost 
because the system is (in general) below SVP. When the 
internal radius of the shell film goes to zero, i.e. y — > 0, 
the difference of grand potentials given by Eq. (3.6) van- 
ishes. The threshold value for the transition SF — ► CC 



is 



Z?*(SF — > CC) = 



g{y} 



l + - 3 Ay 



(3.7) 



provided that y is taken at the minimum of f2g F . 

Before describing the general phase diagram of the sys- 
tems for any A > 0, it is useful to examine with some de- 
tail data calculated for a fixed value of A. Figure 2 shows 
results obtained for A = 0.2 by varying D* and R*. In 
this case some features of and Og F are explicitly in- 
dicated. Solid lines separate domains of stable phases. 
Note that there is a "triple point" below which the space 
is empty (E), to the upper right of which there is a SF 
region, and to the upper left there is CC. The dashed 
straight line in the SF region is the threshold given by 
Eq. (3.4), while dashed curves in the CC regime denote 
the limit for SF solutions determined by the threshold 
given by Eq. (3.5) and the disappearance of the mini- 
mum of Og F . So, the hatched zones indicate parameter 
regions where both CC and SF are negative. However, 
only that phase with lower f2* is stable, the other one is 
metastable and plays an important role in the cycle of 
hysteresis to be addressed later in the paper. 

Let us now turn to the general phase diagram. The re- 
sults for several values of A are displayed in Fig. 1 , where 
only boundaries between stable phases are indicated. In- 
creasing D* favors SF or CC phases against the E phase. 
Which of the condensed phases is stable depends on R*. 
For large R*, SF is typically favored because of the cost 
of CC (the A term) becomes large relative to the benefits 
(from the potential and the decrease of surface energy). 
The SF to E transition curve is rather insensitive to R*. 

Figures 3 and 4 show the reduced phase diagrams pro- 
jected onto the R* - A plane for D* = 5.84 and 11.87 
(i.e. for the 4 He/Cs and 4 He/Na cases, respectively). 
However, for other values of D* the shape of the corre- 
sponding curves is similar. The curves indicate the tran- 
sitions between stable phases. These plots also include 
threshold values obtained from OP-NLDF calculations 
for R p = 10, 12, 15, 20, and 30 A. For 4 He/Na these re- 
sults lie in the neighbor of the "triple point" determined 
by the joining of all three E, SF, and CC phases, where 
the solutions are sensitive to changes of A. Note that the 
E -> CC and SF -> CC transitions provided by the SFM 
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occur at almost the same values of A as that determined 
from OP-NLDF. There is no SF phase for 4 He/Cs sys- 
tems in the analyzed range of R*. On the other hand, 
the predictions for the E — > SF transition do not match 
so well. All these features of the agreement are similar to 
that previously found for planar and cylindrical systems 
(see Fig. 4 in [3] and Fig. 6 in [8]). 



C. Correction to the bare simple fluid model 

The prediction given by the bare SFM for the E — > 
SF phase transition may be improved by introducing a 
correction to the surface term in figF- By looking at 
the surface energy term in Eq. (2.28), which is propor- 
tional to (Ji v , one realizes that it does not vanish in the 
limit of a zero-thickness film (i.e. for £ — > 0). In or- 
der to eliminate this shortcoming one may follow an idea 
adopted by Cheng et al. [34] in writing their Eq. (2.4). 
Accordingly, we shall assume that the surface contribu- 
tion grows exponentially from zero at I = to the bare 
value 4:Trai v [(Rp — r m ) 2 + (R p — r m — £) 2 ] over a charac- 
teristic length ( 

f2 SF (surf) = 4 7rcr /t , [(R p ~ r m f + (R p - r m - tf] 



1 - exp(-£/C) 



This leads to the reduced version 



(3.8 



f2g F (surf) 



1+1 



R* - 1 



1 — exp (— (3 x) 



with 



= r m /C ■ 



(3-9) 



(3.10) 



Let us now analyze the changes of the energetics due 
to this cut-off factor. As expected, its effect only be- 
comes important for rather thin films. In the frame of 
the improved SFM the E — > SF transition is reached for a 
larger value of /io — fx than in the original version. Both 
these features can be observed in Fig. 5, where results 
for 4 He adsorbed into a spherical Na cavity with radius 
R p = 30 A are shown. Hence, the more elaborated ap- 
proach yields a better agreement with OP-NLDF results. 
In fact, for the 4 He/Na systems a very good agreement 
is obtained with (3 a = 3.3 as shown in Fig. 4. It should 
be noted that the corrected curves for the E — > CC and 
SF — > CC transitions cannot be distinguished from the 
previous ones on the scale of the drawing. 

It is plausible to assume that the characteristic length 
£ be mainly proportional to the asymptotic width Woo 
of a free 4 He surface. So, one may write 



where W* is the dimensionless asymptotic width 
W* = W x /r m . 



(3.12) 



By using the asymptotic value Woo — 6 A [33], for 
4 He/Na systems one gets W* ~ 1.3, which leads to 
A a ~ 4.3. A critical film thickness l a may be defined 
by requiring that the exponential of Eq. (3.8) be re- 
duced from unity to about 0.1, i.e. so that the surface 
energy reach a 90% of its bare value of Eq. (2.28). This 
choice for the decay coincides with the density fall off 
adopted for defining the surface thickness in the litera- 
ture [33]. Following this procedure, the obtained value of 
A a yields x a /W* = £ a /Woo — 0.54. This means that, in 
the improved SFM the bare value of the surface energy 
will be essentially reached when an adsorbed film of uni- 
form density po would have approximately the number of 
helium atoms enough to develop a surface with thickness 
Woo at its internal face. 

In Figs. 6 and 7 we compare the chemical potentials, 
the free energies and grand free energies per particle 
provided by the improved SFM with results obtained 
from the OP-NLDF calculations. The data correspond 
to 4 He adsorbed into cavities of Cs and Na with radius 
R p = 30 A. These plots show a good qualitative agree- 
ment between both theoretical approaches, of coarse, 
there are some quantitative differences. Note that in the 
case of 4 He/Na a Maxwell construction is required in or- 
der to determine /j, — /j, = 3.7 K from OP-NLDF data 
for the SF — ► CC transition. 

An example of the density profile evolution as a func- 
tion of the number of 4 He atoms is shown in Fig. 8. 
There OP-NLDF results, including the metastable (or 
unstable) solutions, obtained for a Na cavity with radius 
R p = 30 A are displayed. One may observe that the 
growth of the density profile is continuous and relatively 
smooth throughout the available space. 

It is worthy of notice that the value j3 a — 3.3 fixed 
above lies just in the parameter region determined by 
Cheng et al. [34] for approaching results given by the full 
NLDF theory in the case of planar geometry. Moreover, 
this /3 a is close to f3 = 3 chosen in that paper for the 
purpose of illustration. Therefore it becomes of interest 
to make a connection to the paper of Calbi et al. [4] on 
4 He confined by two parallel walls separated by a dis- 
tance L. By looking at Figs. 15 and 16 of that work one 
may realize that predictions for the E — > F (film) phase 
transitions provided by the bare SFM do not reproduce 
quite well results given by the NLDF theory. So, we ex- 
amined the extend to which a correction similar to that 
included in Eq. (3.8) can improve the agreement in the 
case of adsorption into slits. For such a geometry, the 
surface contribution to the reduced grand free energy is 
given by the factor 4 in Eq. (8) of Ref. [3] . Consequently, 
that term was rewritten as 



(3 = r Tn /( = \r m /W 00 = \/W* 



(3.11) 



Op (surf) = 4 



1-exp - — 



Ax 
W~* 



(3.13) 



G 



here W* is 



(3.14) 



All the reduced dimensionless quantities needed to treat 
a planar system can be obtained from Eqs. (l)-(6) in [4]. 
To illustrate the effect of this cut-off factor we selected 
the case of 4 He adsorbed between planar walls of Li. The 
evaluation was performed with the same cut-off parame- 
ter A a ~ 4.3 utilized for spherical cavities. Figure 9 shows 
data reported in Fig. 15 of [4] together with the shift of 
the E — > F transition towards larger A. From this plot 
one may conclude, that also for the planar geometry the 
agreement with NLDF results is significantly improved. 
As before the E — > CC and SF — > CC transitions remain 
almost unchanged. 



D. Critical radii and adsorption process 

Let us now refer to the phenomenon of hysteresis inher- 
ently related to capillary condensation. For the analysis 
to be made in this section it is convenient to define a 
reduced inner radius r\ as 



V 



y 



R* - 1 



1 



R* - 1 



(3.15) 



Since the empty fraction of the cavity is rf 1 for r\ = 
capillary condensation with a completely filled pore oc- 
curs, while in the r\ — > 1 limit for large R* there are 
very narrow shell films adsorbed on the substrate wall. 
We determined two families of critical inner radii: r\u at 
which the solution for SF becomes metastable and rj C at 
which SF becomes unstable. Associated to the critical 
radii tjm and rjc there are critical reduced chemical po- 
tentials Am and Ac- Numerous experimental evidences 
reported in the literature indicate that porous materials 
fill and drain at different values of the chemical potential 
leading to loops of hysteresis (see Lilly and Hallock [35] 
and references quoted therein). In order to explain these 
loops it was suggested by Cole and Saam [10] that the 
process of adsorption and desorption of a single cavity do 
not always follow the path guided by stable phases. The 
idea is that physical systems exhibit a sort of memory 
trying to remain in the initial phase when A is changed. 

The filling and draining of a spherical cavity may be 
followed by looking at Fig. 6 of Ref . [35] . This process 
can be analyzed by examining Figs. 4 and 10. For a fixed 
pore size and a given adsorbate-substrate combination, 
the adsorption starts at A > in an E phase, which 
corresponds to the vapor phase of the fluid. Next, the 
filling of the pore by decreasing A increases the vapor 
density (this behavior is not included in the SFM, in 
which all the gas states are considered E phase) . At some 
value of A the SF phase becomes stable and a thin shell 
film is formed, then it grows increasing its thickness x and 
reducing the inner radius r\ up to the crossing of fig F with 



hysteretic loop [10] is that after this crossing the filling 
follows the metastable SF phase up to the critical point 
(Ac, xc) where it becomes unstable. Then the SF — > 
CC transition must occur as indicated in Fig. 8(b) of [8]. 
By draining the transition CC — ► SF occurs at (Am, t)m) 
giving rise to an elementary loop of hysteresis. If R* 
is sufficiently small the SF phase cannot exist and the 
filling jumps directly from E to CC following a simpler 
path and no hysteresis takes place. In fact, the whole 
process described here is equal to that assumed in the 
case of cylindrical pores [8,10,35]. 

Figure 10 shows results obtained for adsorption of no- 
ble gases into cavities of Na and Li, for these adsorbate- 
substrate systems D* is always larger than 2.9. For 
each system the nearly universal behavior starts at about 
R* — 1 ~ 15, which for D* in the range 4 — > 20 yields for 
the abscissa 



R* - 1 



(3.16) 



values between 7.5 and 3.6. The limiting curves for 
metastable and unstable SF solutions may be derived 
from expansions in powers of l/V D*. Upon keeping the 
leading term only, the explicit formula for the metastable 
border reads 



* M 



6r? 4 + 1 g 2_g 3 

+ — m 



2\3 



(r 


+ v\ 







while for the unstable one becomes 

18 n 3 



c 



(i-f] 2 y 



(3.17) 



(3.18) 



n 



cc- 



The key assumption for building an elementary 



These functions are also plotted in Fig. 10, where one 
may observe that its behavior is similar to that found 
before for cylindrical pores in terms of hypergeomctric 
functions (see Fig. 10 in [8] and Fig. 2 in [10]). The 
limiting curve corresponding to E — > SF lies at r\ = 1. 



IV. FINAL REMARKS 

A SFM is used for studying adsorption into isolated 
spherical cavities. This model has previously been suc- 
cessfully applied for analyzing phase transitions in the 
cases of planar slits [3] and cylindrical pores [8]. For 
the spherical geometry treated in the present work the 
SFM yields a universal description of capillary conden- 
sation transitions for noble gases confined in pores of al- 
kali metals providing an interpretation of the wide range 
of behavior which can occur. This universality is due to 
the fact that the integrated adsorption potential <?{£} ex- 
pressed in dimensionless quantities [see Eqs. (2.22) and 
(2.26)] is a universal function of £ = r/r m - ln . This prop- 
erty is a consequence of the continuous substrate approx- 
imation. Planar slits and cylindrical pores also exhibit 
this feature, but in the last geometry to get g{£} one 
must perform numerical integrations [8]. 
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The phase diagram for spherical cavities displayed in 
Fig. 1 is qualitatively equivalent to that previously ob- 
tained in for planar and cylindrical geometries (see Fig. 
2 in [3] and Fig. 3 in [8]).' 

The reliability of the SFM was checked by performing a 
comparison with results provided by the much more real- 
istic OP-NLDF. In the latter case the calculations were 
carried out for increasing number of 4 He atoms within 
spherical cavities of different radii. It was found that for 
the E — > CC and SF — > CC phase transitions there is 
a rather good agreement between both approaches (see 
Figs. 3 and 4), while for the E — > SF one there is a size- 
able difference as indicated by Fig. 4. 

It is shown that the bare SFM for shell films given by 
Eq. (2.28) may be improved by including a correction to 
the surface contribution. As can be seen in Figs. 4, 5, 
and 7, the effect of this improvement becomes important 
for thin films. These plots show that the transitions be- 
tween stable phases occur at virtually the same value of 
A in the OP-NLDF theory and in the improved SFM. 
Furthermore, the application of the latter approach for 
planar slits also yields better agreement with NLDF re- 
sults as shown in Fig. 9. 

According to our calculations almost all the examined 
adsorbate-substrate combinations will form CC at T = 
T t , the exceptions are Ne, Ar, Kr, and Xe adsorbed into 
cavities of Cs. In these cases it becomes necessary an 
extra determination of critical temperatures for CC like 
it has been done in [8]. In this respect, it should be 
mentioned that it has been previously found that neither 
Ne confined by planar walls of Cs [3] nor Ne and Ar 
adsorbed into cylinders of Cs [8] form CC at T = T t , 
showing a systematic dependence on the pore curvature. 

The critical radii for metastable phases and unstable 
films displayed in Fig. 10 resemble the behavior found 
previously for cylindrical pores [8,10]. So, adsorption 
potentials built up for different curved geometries lead 
to qualitatively similar features. On the other hand, it 
is suggested how an elementary hysteretic cycle could 
be constructed by filling and draining the cavity along 
paths determined by the critical points (Am,Wm) an d 
(A c , rye). 

In summary, we can state that the present work close 
successfully the application of the so called SFM for 
studying CC in pores exhibiting standard regular geome- 
tries, i.e. planar slits, cylindrical pores and spherical cav- 
ities. 
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TABLE I. Experimental values of relevant observables for 
the inert gases in the liquid phase at the triple point and the 
lattice parameters of solid alkali metals. 



System 


Tt [K] 


Mo [K] 


po [A- 3 ] 


c-iv [K/A 2 ] 


Ref. 


xie 


n a 


















0.274 ± 0.003 


[24,25] 










0.257 ±0.001 


[26] 










0.272 ± 0.002 


[27] 


Ne 


24.55 


-232. 


0.03694 


3.98 


[1,28,29] 


Ar 


83.81 


-930. 


0.02117 


9.74 


[1,28,29] 


Kr 


115.76 


-1342. 


0.01785 


11.22 


[1,28,29] 


Xe 


161.39 


-1907. 


0.01411 


12.65 


[1,28,29] 






a [A]" 


Psub [A" 3 ] C 




Ref. 


Cs 




6.045 


0.009054 




[22] 


Rb 




5.585 


0.0115 




[22] 


K 




5.225 


0.0140 




[22] 


Na 




4.225 


0.02652 




[22] 


Li 




3.491 


0.04701 




[22] 



a The data for 4 He correspond to T = K. 
6 From Ref. [22]. 

c Body centered cubic crystal structure. 



TABLE II. Values of the LJ parameters for the interac- 
tion between noble gases and alkali metals, together with the 
asymptotic well depth and the location of the minimum cor- 
responding to potentials given by Eq. (2.5). 



System 


SLJ [K] a 


0~LJ [A]° 


r m [A] 


£>oo [K] 


D* 


He-Cs 


1.21 


6.47 


5.55 


6.55 


5.84 


He-Na 


1.73 


5.40 


4.63 


15.95 


11.87 


He-Li 


1.92 


5.22 


4.48 


28.34 


20.39 


Ne-Cs 


8.65 


5.23 


4.49 


24.73 


2.06 


Ne-Na 


11.94 


4.37 


3.75 


58.34 


4.06 


Ne-Li 


13.33 


4.23 


3.63 


104.7 


7.06 


Ar-Cs 


51.79 


4.86 


4.17 


118.8 


2.15 


Ar-Na 


60.32 


4.20 


3.60 


261.6 


4.10 


Ar-Li 


66.00 


4.08 


3.50 


465.2 


7.08 


Kr-Cs 


87.16 


4.75 


4.08 


186.7 


2.42 


Kr-Na 


94.10 


4.16 


3.57 


396.6 


4.51 


Kr-Li 


100.10 


4.06 


3.48 


695.2 


7.71 


Xe-Cs 


117.2 


4.87 


4.18 


270.6 


2.52 


Xe-Na 


126.0 


4.27 


3.66 


574.3 


4.70 


Xe-Li 


135.8 


4.16 


3.57 


1015 


8.08 



a Parameters calculated with data taken from [20] by 
applying the procedure outlined in [21]. 
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FIG. 1. Phase diagram showing the behavior at various 
degrees of undersaturation, expressed in terms of the reduced 
chemical potential difference A given by Eq. (2.24). For the 
case A = 0, the line satisfies Eq. (3.1); all values above the 
curve correspond to capillary condensation, while those below 
are "empty". For the other cases, A > 0, there is a "triple 
point" below which the space is empty, to the upper right 
of which there is a shell film, and to the upper left there is 
capillary condensation. 
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FIG. 2. Universal "phase diagram" showing regimes of 
empty, capillary condensed, and adsorbed shell film as a func- 
tion of reduced pore radius and well depth defined in the text 
[see Eqs. (2.11) and (2.21)]. The displayed curves correspond 
to A = 0.2. The hatched zones are regions where both Og F 
and £Iq C are negative. 
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FIG. 3. Reduced phase diagram for D* = 5.84 (i.e. 4 He/Cs 
in our approach). The curve indicates the prediction of the 
SFM for the boundary between possible stable phases (empty 
and capillary condensation) . The full circles stand for transi- 
tions determined from OP-NLDF calculations. 
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FIG. 4. Reduced phase diagram for D* — 11.87 (i.e. 
4 He/Na in our approach). The solid curves indicate the 
prediction for the boundaries among possible stable phases 
(empty, shell film, or capillary condensation) according to the 
SFM. Note the existence of a "triple point". The full circles 
stand for transitions determined from OP-NLDF calculations. 
The dashed curve stand for the E — * SF transition given by 
the improved SFM. 
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FIG. 5. Energetics of 4 He adsorbed into a Na cavity with 
R p — 30 A. Upper part: The grand free energy per particle as 
a function of the number of particles. Lower part: Same for 
the energy per particle and the chemical potential. The solid 
curves are results of the original SFM while dashed curves are 
provided by the improved version. 



FIG. 7. Same as Fig. 6 for 4 He adsorbed into a Na cavity. 
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FIG. 8. Density profiles of 4 He adsorbed into a cavity of 
Na (Rp = 30 A) for N = 200, 400, 600, 800, 1000, 1200, 1400, 
1600, 1630, 1640, and 1800. Dashed curves are metastable or 
unstable solutions. 



FIG. 6. Energetics of 4 He adsorbed into a Cs cavity with 
R p = 30 A. Upper part: The grand free energy per particle as 
a function of the number of particles; Lower part: Same for 
the energy per particle and the chemical potential. The curves 
indicate results provided by the improved SFM. The symbols 
stand for values obtained from OP-NLDF calculations. 
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FIG. 9. Reduced phase diagram for 4 He confined by two 
planar walls of Li (D* = 12). Solid curves are SFM predic- 
tions and full circles DF results, both of them taken from [4] 
(see text). The dashed curve shows how much the correction 
introduced in the present work improves the prediction of the 
SFM for the E — > F phase transition. 



FIG. 10. Nearly universal correspondence of critical values 
of reduced inner radii r) M and 7]c and the dimensionless ratio 
of effective radius with potential strength (R* — D* in 
the single hysteretic loop. The solid curves are asymptotic 
functions given by Eqs. (3.17) and (3.18). 
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